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ABSTRACT 

This  report  considers  empirical  solutions  to  the  first-passage 
or  single-highest-peak  (SHP)  problem  for  a  distributed  elastic  structure 
with  rectangular  geometry  subjected  to  both  stationary  and  a  specific 
form  of  nonstationary  random  excitation.  The  structure  is  a  flat,  homo¬ 
geneous,  uniform,  square  plate  and  the  applied  stationary  excitation 
is  white  noise  perfectly  correlated  in  both  space  and  time.  The  non¬ 
stationary  excitation  is  a  rectangular  noise  burst  with  a  unity  correla¬ 
tion  in  both  space  and  time.  The  structure  and  the  excitation  are  simu¬ 
lated  electrically  and  peak  response  data  are  collected  for  (1)  simply 
supported  boundary  conditions  at  all  edges  and  (2)  rigidly  clamped 
boundary  conditions  at  all  edges.  These  response  data  are  used  to 
establish  probability  curves  yielding  an  estimate  of  the  probability  that 
the  maximum  response,  for  a  finite  time  interval,  remains  below  a  pre¬ 
selected  threshold  level. 
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1.  INTRODUCTION 


The  SHP  (single-highest-peak)  problem  generally  refers  to  predict¬ 
ing  the  maximum  response  a  mechanical  system  may  achieve,  within  a 
finite  time  interval,  when  subjected  to  stationary,  broad-band  white  noise 
excitation.  It  is  noted  that  the  maximum  response  need  not  necessarily 
correspond  to  a  peak  amplitude  in  the  response  time  history.  The  solu¬ 
tion  rests  with  determining  P-^(|P|  <  Pq  »  T)  which  is  the  probability  that 
the  absolute  value  of  the  response  p  remains  below  a  threshold  level  p 

o 

within  the  time  interval  T  .  The  term  p  is  a  dimensionless  ratio  of  the 
maximum  amplitude  to  RMS  response  of  the  system  and  T  is  the  time 
interval  over  which  the  response  is  observed.  The  p  response  maybe 
a  positive  maximum  value  p+,  a  negative  maximum  value  p  ,  or  an 
absolute  maximum  value  IM  and  the  RMS  response  is  the  system  RMS 
response  to  stationary  white  noise  excitation  perfectly  correlated  in 
space  and  time.  For  this  discussion,  only  solutions  involving  |p| 
values  are  shown. 

For  a  single  degree-of-freedom  system,  the  SHP  problem  has 
been  examined  in  detail  by  various  investigators  using  both  simulation 
and  analytical  procedures  as  is  discussed  in  References  listed  at  the 
end  of  this  report.  Solutions  acceptable  for  practical  applications  have 
been  cleverly  determined  not  only  for  stationary  white  noise  excitation 
(References  1  and  7),  but  also  for  white  noise  applied  as  a  step  function 
(Reference  5)  and  for  particular  noise  bursts  as  well  (References  4  and  8). 
It  is  judicious  to  examine  some  of  the  established  results  for  the  vibra¬ 
tion  of  a  single  degree-of-freedom  system  before  discussing  the  SHP 
problem  for  a  distributed  elastic  structure. 
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2.  RESULTS  FOR  A  SINGLE  DEGREE-OF-FREEDOM  SYSTEM 


A  linear  single  degree-of-fr  eedom  mechanical  system  (also  called 
a  linear  mechanical  oscillator)  is  shown  as  Figure  1 


I 

f(t)  I 


Figure  1.  A  Single  Degree-of-Freedom  Mechanical  System 

The  dynamic  behavior  of  this  system  is  defined  by  the  equation  of  motion 

mx  +  c'x  +  kx  =  f(t)  (1) 

where  m  corresponds  to  the  mass,  c  the  viscous  damping  coefficient, 
k  the  spring  constant  and  f(t)  a  forcing  function  applied  to  the  mass. 

In  an  alternate  form,  this  equation  may  be  written  as 

x  +  2  £  go  x  +  go  x  =  (2) 

n  n  m 
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where 


—  =  2£co 

m  n 

(3) 

—  =  w2  =  (2trf  )2 

m  n  n 


and  £,  is  the  damping  factor  and  w  the  undamped  natural  frequency  of 

the  system.  For  f(t)  equal  to  a  simple  harmonic  forcing  function, 

say  f(t)  =  f  sin  wt,  the  forced  displacement  response  x(t)  is 
o 

x(t)  =  — f  sin  cot  (4) 

c \  o 

moo 

n 

where  the  magnification  factor  H(co)  is  defined  in  complex  form  by 


H(co)  = 


1 


1 


2 

+  i  2£  — 
co 
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(5) 


For  f(t)  equal  to  an  arbitrary  deterministic  function  of  time,  the  re- 
sponse  may  be  calculated  in  the  following  ways 


H(to)  F(w)  elwt  dto 


h(T)  f(t  -  t)  dr 


(6) 
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where  F(co)  is  the  Fourier  Transform  of  f(t)  and  h(r)  is  the  system 
response  to  a  unit  impulse  excitation.  The  solution  obtained  by  inte¬ 
grating  over  oo  is  a  Fourier  integral  solution  while  the  form  given  by 
integrating  over  t  is  a  convolution  integral  solution.  The  Fourier 
integral  gives  the  solution  x(t)  as  a  superposition  of  steady-state  re¬ 
sponses  and  the  convolution  integral  gives  the  solution  as  a  super¬ 
position  of  free  vibration  responses. 

The  output  response  of  a  lightly  damped  single  degree- of-freedom 
system  wherein  f(t)  is  stationary  broadband  white  noise  excitation  is 
shown  as  Figure  2.  The  undamped  natural  frequency  of  the  mechanical 
system  is  assumed  large  compared  with  the  half-power  bandwidth  AND 
the  bandwidth  of  the  input  spectrum  is  assumed  wide  in  comparison  to 
the  system  half-power  bandwidth  and  includes  f  .  It  is  further  implied 
the  system  has  achieved  stationarity  in  its  response  and  we  are  ob¬ 
serving  the  response  beginning  at  an  arbitrary  instant  of  time. 


Amplitude 


Time 


Figure  2.  Response  of  a  Lightly  Damped  Mechanical  Oscillator  to 
Broadband  Stationary  Random  Excitation 
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This  response  appears  as  a  sinusoid  (at  the  frequency  f  )  with 

n 

a  slowly  varying  random  amplitude  and  random  phase.  The  average 
number  of  zero  crossings  per  second  E  [N  ]  is  approximately 


E[N  ]  =-' 

O  IT 


2  2 
W1  +0)lw2  +  W2 


00 


n 

IT 


(?) 


where  oo^  and  oo^  are  interpreted  as  the  half-power  point  frequencies 
in  radians  per  second.  The  mean  square  displacement  response  ^  is 


7  ttG  Q 
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x  Z  5  n 

2  m  oo 
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where 


Q  = 


2£ 


and  is  the  magnitude  of  the  input  spectrum  with  units  of  lb  /rad/ sec. 
The  term  I  (see  page  80  of  Reference  6  and  Reference  11)  is  a  dimension¬ 
less  coefficient  ranging  from  zero  to  one  with  its  exact  value  being  de¬ 
pendent  upon  (1)  the  damping  factor  £  and  (2)  the  ratio  of  the  cutoff 

frequency  oo  of  the  input  spectrum  to  oo  as  shown  in  Figure  3.  Simi- 

c  n  2 
larly,  the  mean  square  velocity  response  i(j%  is 
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Figure  3.  Numerical  Value  of  I 


where  the  dimensionless  coefficient  II  is  shown  as  Figure  4.  For 

high  Q  systems  (L  <  0.  010)  ,  both  I  and  II  are  nearly  zero  for  go  /go 

n  n  c  n 

somewhat  less  than  one;  they  rapidly  approach  unity  as  go  /go  becomes 
slightly  greater  than  one  and  converge  to  unity  as  go  /go  — oo  .  This  be¬ 
havior  implies  the  mechanical  system  acts  as  a  highly  selective  band¬ 
pass  filter  for  both  displacement  and  velocity  and  admits  frequency  com¬ 
ponents  principally  near  go  .  By  selecting  I  and  II  as  unity,  upper 
*  2  2  n  n  n 

limits  for  ib  and  ib »  are  calculated, 
x  x 

The  autocorrelation  function  for  the  displacement  response  is 
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where  go^  is  the  damped  natural  frequency  of  the  system  and  related 
to  go^  as  go^  =  GOnYT!  .  The  magnitude  of  this  function  at  t  =  0  yields 
the  mean  square  displacement  response 
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which  is  identical  to  Eq.  (8)  with  1^  equal  to  unity. 


❖ 

The  mean  square  acceleration  response  may  be  developed  in  a  similar 

manner  and,  in  general,  tends  to  be  divergent  as  go  /go  >1  and  becomes 

unbounded  as  to  /go  —*oo.  1  11 

c  n 
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Figure  4.  Numerical  Value  of  II 


The  time  lapse  At  for  the  amplitude  of  R^(t)  to  decay  to  —  times  its 
initial  amplitude  is  given  by 


At 


Q 

irf 

n 


1 

n 


(ID 


and  is  called  the  correlation  time  or  correlation  interval  (Reference  5). 
This  time  interval  corresponds  to  the  time  lapse  for  the  free  vibration 
of  a  lightly  damped  single  degree- of-freedom  system  to  decay  to  -■ 
times  an  initial  amplitude  . 

For  a  white  noise  excitation,  the  probability  density  function  of 
the  response  at  any  instant  of  time  is  Gaussian  (Normal).  The  probability 
density  of  the  envelope  of  the  peaks  is  Rayleigh.  The  probability  density 
of  the  individual  peak  amplitudes  is,  likewise,  nearly  Rayleigh  and  be¬ 
comes  precisely  Rayleigh  when  the  damping  becomes  sufficiently  small 
to  define  an  infinitely  narrow  bandpass  at  f  , 

The  probability  density  of  the  response  maxima  during  a  finite  ob¬ 
servation  time  is  neither  Gaussian  nor  Rayleigh  and  forms  the  crux  of  the 
SHP  problem.  A  solution  for  this  probability  distribution  function  is 
given  in  Reference  10  as 


M 


<IM  <  P„ •  T)  =  A  e 


-a  Q oo  T 
o  n 


T  >  T 


cor 


(12) 


where  A  is  a  constant  whose  magnitude  is  dependent  upon  the  initial  con- 
o 

ditions  of  the  system,  and  t  denotes  the  time  lapse  for  the  autocorre- 
3  cor 

lation  function  to  decay  to  a  negligible  value.  For  most  applications  with 
(3^  >  2  ,  A^  may  be  assumed  equal  to  unity  with  trivial  error.  The 
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EXPONENTIAL  CONSTANT 


term  a  is  dependent  upon  the  system  damping  Q  and  a  preselected 
o 

level  of  the  6  ratio.  Plots  of  a  versus  Q  are  shown  in  Figure  5 
°  o  . 

as  families  of  curves  in  (3^  for  both  (3  and  |p|  .  This  figure  in  con¬ 
junction  with  Eq.  (12)  provides  a  theoretical  solution  to  the  SHP  prob¬ 
lem  for  a  linear  mechanical  oscillator. 

Empirical  solutions  to  this  problem  have  been  determined  using 
techniques  amenable  to  digital  and  analog  computation.  By  means  of  an 
active  electrical  analog  simulation  (differential  amplifier  circuits),  the 
results  shown  as  Figures  6,  7,  8  and  9  have  been  obtained.  Figures  6 
and  7  are  families  of  plots  in  Q  of  the  dimensionless  time  parameter 
f  T/Q  versus  |(3|  for  PM(|P|  <  PQ)  =  0.50  and  0.95.  These  data  may 
be  compared  directly  with  the  results  of  Figure  5  by  restating  Eq.  (12)  in 
the  form 


i»[fm<|p|  <  e0>] 


(13) 


and  evaluating  this  expression  in  the  following  manner 


These  data  are  partially  from  Reference  4  and  partially  from  inde¬ 
pendent  analog  studies  carried  out  after  the  writing  of  Reference  4.  Addi¬ 
tional  data  were  obtained  from  a  digital  simulation  study  sponsored  by 
Mr,  E.  H.  Schell  of  AFFDL,  Wright- Patters  on  Air  Force  Base,  Ohio. 
These  digital  results  tend  to  substantiate  the  shown  analog  data  although 
the  (3  results  for  the  higher  probability  values  were  consistently  lower 
(within  10%  and  generally  not  more  than  5%). 
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The  data  values  for  6  =3  and  6  =  4  in  Figure  6  and  for  6=3, 

o  o  °  o 

=  4  and  (3^  =  5  in  Figure  7  agree  closely  with  their  equivalent  values 
in  Figure  5.  It  may  thus  be  concluded  that  these  analog  results  defini¬ 
tively  support  the  theoretical  solution  of  Eq.  (12). 

Figure  8  is  a  set  of  plots  in  Q  of  |p|  versus  f  T/Q  for 
P^(|  P|  <  Pq)  =  0.  95  and  P^(|  p|  <  p  )  =  0.  50.  Only  curves  for  Q  =  20 
and  Q  =  50  are  shown  since  these  values  are  considered  to  be  more 
typical  of  damping  in  physical  structures,  and  thus,  of  most  interest. 
These  curves  are  exploded  views  of  Figures  6  and  7  for  f  T/Q  <  5  and 
find  application  wherein  the  time  interval  T  may  be  relatively  short. 

The  larger  |p|  values  are  associated  with  the  lower  damping  ratios  for 
all  fRT / Q  with  the  sharpest  rate  of  increase  noted  for  f  T/Q  ranging 

from  0  to  3  .  It  is  observed  that  P  values  of  3  for  P  (Ipl  <  p  )  =  0.  50 

^  o  M'  '  —  ro 

and  p^  values  of  4  for  P^(|  p|  _<  p^)  =  0.  95  are  to  be  expected  where 

f  T/Q  >  5. 
n  — 

Figure  9  shows  |p|  versus  f^T/Q  plots  of  the  average  maxima  of 
| P |  for  Q  =  20  and  Q  =  50  .  These  curves  are  very  nearly  the  same  as 
those  for  P^.(|p|  <  P^)  =  0.  50  in  Figure  8  with  one  exception.  The  data 
for  Q  =  20  appears  to  become  numerically  equivalent  with  that  for  Q  =  50 
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Figure  8.  Maximum  Response  of  a  Single  Degree -of-Freedom  System  to  Stationary 
White  Noise 
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to  Stationary  White  Noise 


where  f  T/Q  >  4.  From  these  data,  the  median  and  average  values  of 
n  — 

|p|  for  Q  =  50  tend  to  become  very  nearly  the  same.  It  is  suggested 
that  a  practical  estimate  of  the  average  maximum  response  may  be  de¬ 
termined  by  using  the  | (3|  versus  f^T/Q  plot  for  Q  =  50  . 

The  maximum  response  of  a  simple  oscillator  to  shaped  noise 
bursts  also  has  been  examined  by  an  active  analog  simulation  as  dis¬ 
cussed  in  References  4  and  8.  The  noise  burst  excitation  may  be  written 
as 

f(t)  =  E(t)  n(t)  (14) 


where  E(t)  is  a  specifically  defined  envelope  function  and  n(t)  is  the 
output  response  of  a  random  noise  generator  with  a  white  noise  spectrum 
of  magnitude  .  Data  from  these  studies  are  shown  as  Figures  10 
and  11  where  IpI  is  the  dimensionless  ratio  of  the  maximum  response 
to  the  RMS  response  and  f  t/Q  is  the  dimensionless  time  ratio  where  t 
is  a  time  duration  of  the  input  burst  of  noise.  The  RMS  response  is  the 
output  response  of  the  oscillator  to  broadband  white  noise  wherein 
E(t)  =  1  .  By  normalizing  the  p  response  in  this  way,  time  varying 
RMS  considerations  are  avoided.  These  data  suggest  the  noise  burst 
response  maxima  become  relatively  the  same  as  that  for  the  stationary 
response  when  f^r/Q  -1.5.  Thus,  plots  such  as  these  may  be  used  to 
establish  criteria  for  the  conditions  under  which  an  oscillator  achieves 
stationarity  in  its  response. 
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Figure  10.  Response  Maxima  of  a  Single  Degree-of-Freedom  System  to  a  Rectangular 
Burst  of  White  Noise 


Response  Maxima  of  a  Single  Degree-of-Freedom  System  to  a  Rectangula 
Burst  of  White  Noise 


3.  RESPONSE  MAXIMA  OF  DISTRIBUTED  ELASTIC  PLATES 


In  this  section,  the  SHP  problem  is  considered  wherein  the  struc¬ 
ture  is  a  flat  square  plate  and  the  applied  force  excitation  is  stationary 
white  noise  correlated  perfectly  in  both  space  and  time.  The  (3  re¬ 
sponse  statistics  are  then  examined  via  measurements  taken  at  a  point 
near  the  center  of  the  plate  for  two  types  of  boundary  conditions,  (1) 
simply  supported  at  all  edges  and  (2)  rigidly  clamped  at  all  edges. 

The  equation  of  motion  for  a  flat  square  plate  is 


+  mw  =  f (x ,  y,  t) 


(15) 


where  f(x,  y,  t)  is  the  applied  force  excitation  per  unit  area,  m  the  mass 
per  unit  area,  w  =  w(x,  y,  t)  =  the  lateral  deflection  from  static  equilibrium, 
v  the  Poisson’s  ratio  and 


r> 


d2(  ) 

dt2 


D  = 


Eh 


12(1  -  v  ) 


(16) 


The  coefficient  D  is  the  flexural  rigidity  of  the  plate  wherein  E  is 
Young's  modulus  and  h  the  plate  thickness.  The  notation  ) 
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defines  a  spatial  operator  suitable  to  the  rectangular  geometry  of  the 
square  plate  where  ^  is  the  familiar  dell  operator  for  rectangular  co- 
ordinates. 

The  mean  square  displacement  response  of  an  arbitrary  linear 
elastic  structure  to  a  uniform,  homogeneous  forcing  field  may  be 
written  as 

/oo 

H  (co)  H*(co)  L  (co)  dco  (17) 

-oo  J  J 


where  the  quantity  L^(co)  is  expressed  in  terms  of  the  spatial  cross- 
spectral  density  function  S  (r,  r',co)  as 


A  -  A 


L  (co)  =  __  — -1  r  [  [  (j).(r)  <t>(r')  S  (?■,“',  co)  dr  dr' 

jk  M.^u.  «kJo  JO  J 


(18) 


or  in  terms  of  the  joint  acceptance  j  (go)  as 

JK 


S  (r  ,  u>)  A 

V"’  =  ^  -  2"  "2  jjk(u) 

M.  M.  co.  oo 
j  k  j  k 


(19) 


The  notation  T  refers  to  a  spatial  position  vector,  $  (F)  the  jth  normal 
mode  of  the  system,  co^  the  jth  modal  frequency,  H^(co)  the  jth  modal 
magnification  factor  defined  as 
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H.(co)  = 
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and  H^co)  the  complex  conjugate  of  H^(oo)  .  The  symbol  A  denotes  the 

area  of  the  plate  and  the  quantity  S^(r^,  go)  refers  to  a  spectral  density 

function  of  the  applied  force  at  the  spatial  location  F  .  If  the  position 

r  is  selected  to  be  a  point  on  the  structure  where  S„(r*  ,  go)  is  a 
°  2  f  o 

maximum,  j  (co)  will  vary  from  zero  to  one.  From  Eqs.  (18)  and  (19), 

jk 

the  joint  acceptance  is  found  to  be 
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By  assuming  the  j  4  k  terms  to  be  negligible  in  comparison  to  the 
j  =  k  terms,  the  mean  square  displacement  response  becomes  (for  the 
conditions  of  this  problem) 


or,  in  terms  of  the  physically  realizable  spectrum  G  , 
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Similarly,  the  mean  square  velocity  response  may  be  written  as 
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In  terms  of  the  functions  I  and  II  shown  as  Figure  3  and  4,  Eqs.  (23) 

n  n 

and  (24)  are  of  the  form 
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The  SHP  problem  for  the  plate  structure  is  studied  empirically 
by  examining  the  (3  results  from  the  electrical  analog  simulation  de¬ 
picted  schematically  by  Figure  12.  The  transformer  whiffle-tree  cir¬ 
cuit  is  used  to  create  a  spatial  correlation  of  unity  for  the  white  noise 
excitation  which  is  applied  to  the  plate  analog.  As  shown  in  this  figure, 
the  data  of  interest  include  both  positive  and  negative  response  maxima 
(measured  by  the  "peak"  readers)  and  the  RMS  response. 
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Excitation  Velocity 

Respons  e 

Figure  12.  Block  Diagram  of  Analog  Simulation  Study 


The  plate  analog  circuits  are  derived  from  energy  concepts  (Ref¬ 
erence  2)  and  may  be  categorized  as  mobility  analogs.  Such  circuits 

consist  of  passive  electrical  components,  are  appropriate  for  nonuniform 

* 

physical  and  geometric  properties,  and  appear  topologically  similar 
to  the  physical  system.  Mathematically,  this  type  of  analog  corresponds 
to  a  finite -difference  model  whereas,  mechanically,  it  may  be  interpreted 
as  a  form  of  lumped-parameter  model.  Background  information  about 
such  structural  analogs  is  found  in  Reference  9- 

The  physical  dimensions  and  physical  properties  of  the  uniform, 
homogeneous  plate  are  defined  as 


For  the  more  complicated  structural  configurations,  a  rather  vivid 
imagination  admittedly  is  required. 
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a  =  b  =  j0  =  llin 


h  =  0.  0625  in 
E  =  10.  5  x  106  lb/in2 
p  =  0.  10  lb /in2 
v  =  0.  3 

For  simply  supported  edge  conditions,  the  plate  modal  frequencies  are 
given  by 


where  rs  denote  the  mode  numbers  (the  number  of  half- sines  J  in  the 

-  5  2  3 

x  and  y  dimensions).  With  D  =  235  lb-in  and  m  =  1.  62  x  10  lb-sec  /in  , 
the  fundamental  modal  frequency  is  ^  99  cps  .  For  rigidly  clamped 
edge  conditions,  the  fundamental  modal  frequency  is  f^  —  181  cps. 

The  plate  circuits  and  expressions  for  the  component  values  used 
in  this  study  are  given  in  Reference  3.  By  scaling  the  analog  circuits 
such  that  electrical  time  corresponds  to  real  time  (N  =  1)  and  applying  a  uni¬ 
formly  distributed  sinusoidal  load,  the  velocity  to  force  frequency  response 
41  41 

functions  at  x  =  — and  y  =  are  obtained  as  shown  in  Figure  13. 

These  response  functions  indicate  the  magnitude  of  the  circuit  impedance 
at  a  particular  point  and  are  equivalent  to  mobility  plots  (magnitude  only). 

The  fundamental  modal  frequencies  are  noted  to  be  approximately  equal 
to  the  values  previously  calculated  for  both  simply  supported  and  rigidly 
clamped  boundary  conditions.  Due  to  the  symmetry  of  the  loading,  only 
even  valued  modes  (i.  e.  where  the  mode  shapes  are  symmetric  about  the 
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VELOCITY  RESPONSE  IN  db  VELOCITY  RESPONSE  IN  db 


FREQUENCY  IN  CYCLES  PER  SECOND 


Figure  13.  Frequency  Response  Functions  of  a  Square  Plate 
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mid-spans  of  the  plate)  can  be  excited  as  is  seen  in  these  plots.  The 
depth  of  the  valleys  between  the  modal  frequencies  is  indicative  of  the 
modal  coupling  in  the  system. 

A  measure  of  the  damping  in  the  first  mode  is  given  by 


Q 


11 


(28) 


where  Af^  is  the  half-power  bandwidth  centered  at  f^  .  This  expression 
yields  Q  —  33.  9  for  the  simply  supported  boundaries  and  —  39.  7 
for  the  rigidly  clamped  boundaries.  Such  damping  is  inherent  in  the  com¬ 
puter  (due  to  the  fact  that  the  circuit  capacitors,  inductors  and  trans¬ 
formers  are  not  ideal,  interconnecting  wires  have  resistance,  physical 
layout  of  the  components  produce  parasitics,  etc.  )  and,  although  it  can¬ 


not  be  precisely  controlled,  damping  values  in  the  first  mode  ranging  from 
Q  =  30  to  Q  =  50  are  to  be  expected  with  proper  scaling. 

In  the  simulation  study,  the  plate  is  continuously  excited  by  station¬ 


ary  white  noise  correlated  perfectly  in  both  space  and  time  and  the  velocity 
4i  4f 

response  at  x  =  —  and  y  =  —  is  observed  for  various  time  durations  T 
9  9 

such  that  the  dimensionless  ratio  f  T/Q^  ranges  from  0  to  5.00. 

The  magnitude  of  the  white  noise  input  is  adjusted  such  that  a  response 
maxima  of  approximately  six  times  the  RMS  response  (6  values)  could 
occur  without  distortion.  For  a  specified  T  value,  the  response  is 
sampled  100  times  wherein  both  the  positive  and  negative  response  max¬ 
ima  are  recorded  for  each  sampling  thus  producing  a  total  of  200  data 
points.  Typical  response  time  histories  (as  observed  over  a  finite  sam¬ 
pling  time  interval  T)  are  shown  as  Figures  14  and  15.  It  is  to  be  noted 


the  response  does  not  appear  to  be  the  same  as  the  stationary  response  of 
a  linear  oscillator  to  white  noise,  i.  e.  ,  as  a  type  of  harmonic  motion  with 
randomly  varying  amplitude  and  phase.  Rather,  higher  frequency  effects 
due  to  modal  frequencies  other  than  f^  are  evident. 
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Boundaries 
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Rigidly  Clamped  Boundaries 


These  individual  maxima  response  data  are  next  arranged  in 
rank  order,  then  divided  by  the  RMS  response  to  form  histograms  of 

-j-  /V 

PjvI(P  <  PQ)  and  PM(P  <  PQ)  •  A  typical  example  of  such  data  is 
shown  as  Figure  16.  By  adding  the  p+  and  p  data,  a  histogram  of 
Pw(lpl  <  Pq)  is  formed  as  shown  in  Figure  17  and  this  is  used  (one 
such  histogram  for  each  sampling  time  T)  to  obtain  individual  prob¬ 
ability  data  points  for  |  (3^  versus  £^  T/Q  plots.  It  is  understood 
that  these  histograms  are  estimates  of  the  true  maxima  response  prob¬ 
ability  density  functions  for  a  specified  value  of  fjjT/Q  .  Additional 
data  should  therefore  be  accumulated  so  that  the  variance  of  such  data 
may  be  analyzed. 

Following  such  an  empirical  procedure  just  described,  the  |p| 
versus  fjjT/Q  plots  of  Figures  18  and  19  are  developed.  Note  that 
only  the  fundamental  modal  frequency  f  and  its  associated  damping 
need  be  known  for  the  structure.  These  curves  display  a  be¬ 
havior  similar  to  that  for  the  single  degree-of-freedom  system  (com¬ 
pare  with  Figures  8  and  9)  but  have  consistently  higher  values  of  | {3 1 

A 

for  a  given  probability.  The  curve  for  (I  p|  <  pc  )  =  0.  50  is  approxi¬ 
mately  the  same  as  that  for  the  average  |p|  and  consequently  is  not  in¬ 
cluded  as  a  separate  figure. 

For  a  given  probability,  the  |(3j^  ratio  increases  indefinitely  with 

time  and  theoretically  | (3 1  — ♦  oo  as  T  oo  ,  the  latter  being  a  truism 

of  limited  practical  concern.  The  greatest  rate  of  increase  in  |p| 

is  seen  to  occur  within  0<f^T/Q^<3.  For  f^T/Q  -  3  ,  |(3| 

values  of  approximately  three  are  expected  on  the  average  and  95%  of 

all  amplitude  maxima  are  expected  to  be  equal  to  or  less  than  (3  ~  4  . 

o 
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Simply  Supported  Boundaries 
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igure  18.  Maximum  Response  of  a  Square  Plate  to  White  Noise  Correlated  Perfectly 
in  Space  and  Time 


Figure  19.  Maximum  Response  of  a  Square  Plate  to  White  Noise  Excitation  Correlated 
Perfectly  in  Space  and  Time 


Consider  the  maximum  response  characteristics  of  a  plate  sub¬ 
jected  to  a  burst  of  white  noise  shaped  in  time  by  a  well  defined  envelope 
function.  It  is  required  to  develop  a  solution  such  that  one  may  predict 
the  expected  single  highest  response  maximum  due  to  a  noise  burst  of 
finite  duration.  This  problem  may  be  examined  in  much  the  same  way 
as  with  the  stationary  excitation,  that  is  by  an  analog  simulation  study 
wherein  data  for  the  response  maxima  are  collected  and  arranged  to 
form  |  p|  plots.  However,  special  attention  must  be  given  to  (1)  the 
RMS  response  used  to  form  the  p  ratio  and  (2)  the  time  duration  com¬ 
mon  to  the  dimensionless  time  parameter  f  T/Q^  .  As  with  the 
previous  problem,  the  p  response  is  normalized  using  the  RMS  response 
of  the  system  to  broadband  white  noise  correlated  perfectly  in  both  space 
and  time  wherein  the  magnitude  of  the  stationary  noise  corresponds  to 
the  value  where  the  envelope  function  is  a  maximum  (E(t)  =  1)  .  Al¬ 
though  admittedly  arbitrary,  this  procedure  avoids  consideration  of  a 
time  varying  RMS.  The  sampling  time  interval  T  is  replaced  by  T 
which  denotes  the  time  duration  of  the  noise  burst.  With  these  modifica¬ 
tions  in  definition,  this  nonstationary  problem  may  now  be  plausibly 
attacked. 

With  the  system  initially  at  rest,  a  rectangular  noise  burst  is 
applied  for  r  seconds  typically  producing  the  response  as  shown  in  Fig¬ 
ure  20.  The  system  response  is  monitored  from  the  onset  of  the  ex¬ 
citation  until  the  free  vibration  response  has  decayed  to  a  negligible  level 
and  both  the  positive  and  negative  maximum  response  values  are  re¬ 
corded.  By  sampling  the  system  response  one  hundred  times  for  each  r 
setting  and  arranging  the  data  in  the  same  manner  as  with  the  stationary 
excitation,  the  |p|  plots  of  Figures  21  and  22  are  constructed.  Data 
points  for  simply  supported  boundaries  are  noted  by  »'s  and  data  points 
for  rigidly  clamped  boundaries  are  shown  by  A' s  . 
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r related  Perfectly  in  Space  and  Time 


Maximum  Response  of  a  Square  Plate  to  Rectangular  Bursts  of  White  Noise 
Correlated  Perfectly  in  Space  and  Time 
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'igure  22.  Maximum  Response  of  a  Square  Plate  to  Rectangular  Bursts  of  White  Noise 
Correlated  Perfectly  in  Space  and  Time 


Similar  to  the  behavior  for  a  single  degree-of-freedom  system, 
the  response  data  are  bounded  by  the  curve  depicting  the  stationary  re¬ 
sponse  of  the  plates.  As  the  time  duration  of  the  noise  burst  becomes 
long  relative  to  the  fundamental  modal  period  of  the  structure,  the  non¬ 
stationary  excitation  appears  as  a  stationary  white  noise  input  and  the 
noise  burst  response  becomes  essentially  that  of  the  stationary  curve. 
Curves  such  as  these  may  be  used  to  define  the  conditions  under  which 
a  multimode  system  excited  by  a  noise  burst  achieves  stationarity  in  its 
response. 

Of  vast  importance,  however,  is  that  the  data  for  both  stationary 
and  nonstationary  random  excitation  imply  the  I  Pi  response  tends  to  be 
independent  of  the  boundary  conditions.  Such  a  statement  must  be  con¬ 
sidered  as  a  qualified  conjecture  until  more  data  are  collected  and  ana¬ 
lyzed  and/or  theoretical  work  is  advanced  to  support  or  reject  this 
claim.  In  the  interim,  |(3|  curves  such  as  included  in  this  report  may 
be  used  discretely  in  engineering  applications  wherein  a  maxima  re¬ 
sponse  criteria  for  distributed  structures  is  applicable. 
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4.  CONCLUDING  REMARKS 


The  SHP  problem  has  been  reviewed  briefly  for  a  single  degree- 
of-freedom  system  wherein  both  analytical  and  empirical  solutions  have 
been  discussed.  More  important,  data  have  been  presented  forming  an 
empirical  solution  to  the  SHP  problem  for  a  square  plate  with  (1)  simply 
supported  boundaries  and  (2)  rigidly  clamped  boundaries.  The  random 
excitation  applied  over  the  plate  is  (1)  stationary  white  noise  with  a 
unity  correlation  in  both  space  and  time. and  (2)  nonstationary  excitation 
shaped  (in  time)  as  rectangular  noise  bursts  with  a  unity  correlation  in 
both  space  and  time.  A  comparison  of  ^^(IPI  i  PQ)  data  for  the  single 
degree-of-freedom  system  and  plate  shows  the  |P|  re  sponse  maxima  to 
be  consistently  higher  for  the  plate.  Of  particular  significance  in  the 
plate  data  is  the  implication  that  the  P^.(|P|  <  P  )  curves  are  inde¬ 
pendent  of  the  boundary  conditions.  Such  an  implication  has  far-reaching 
practical  significance  and  should  be  investigated  thoroughly  in  the  im¬ 
mediate  future.  In  addition,  the  effect  on  the  response  maxima  of  mass 
loading  and  of  bandlimited  excitation  with  spatial  and  spectral  character¬ 
istics  more  typical  of  service  conditions  would  be  indeed  worthwhile 
investigations. 
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